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SUMMARY 


Microcomputer-based  algorithms  for  lognormal  distribution  parameter 
estimation  are  presented.  A  program  diskette  will  be  made  available  upon 
request.  For  given  data  the  parameters  (three  or  two)  are  estimated  (i)  by  means 
of  the  moments,  and  (ii)  by  means  of  the  maximum-likeiihood  principle. 
Maximum-likelihood  estimation  is  done  in  two  different  ways:  (i)  by  means  of 
the  derivative  equations  which  result  from  the  logarithmic  likelihood  function, 
and  (ii)  by  means  of  direct  optimization  of  the  logarithmic  likelihood  function 
itself.  Moment  estimates  are  used  as  starting  values  for  the  optimization 
process. 
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I.  INTRODUCTION 


The  lognormal  probability  distribution  class  is  defined  by  its  probability  density 
function  (PDF) 


f(x)  = 


exp 


[log(x-X)-;i]^  .x>x 


^0  .  X  S  X  . 


(") 


The  distribution  parameters  are^  shift  X  e  R  ,  shape  a  >  0,  and  scale  p  e  R.  The 
scaling  property  of  p  becomes  evident  under  the  transformation  p  =  log  2f  ,  >  0,  so 

that  log  (x-X)  -  p  =  log  [(x-X)2f~M. 

The  lognormal  PDF  (•*)  is  unimodal.  It  takes  exactly  one  maximum  in  the  internal 
(X.oo)  which  is  located  at  x  =  X  +  exp  (p-a^)  (mode).  It  is  of  diffusion  type,  i.e.,  it  is 
generated  by  a  one-dimensional  parabolic  differential  equation  of  Fokker-Planck  type. 

Two  algorithms  for  estimation  of  the  three  parameters  of  (*<)  relative  to  given 
statistical  data  are  presented  which  are  based  on  the  maximum-likelihood  (ML) 
principle.  Both  can  be  used  for  sample  (raw)  data  as  well  as  histogram  (grouped)  data. 
One  of  them  utilizes  the  three  equations  which  are  obtained  by  equating  to  zero  the 
three  first  partial  derivatives  of  the  associated  logarithmic  likelihood  function  (LLi^). 
Two  of  the  parameters  (shape  and  scale)  can  be  eliminated  so  that  only  one  equation, 
h(X)  =  0,  actually  remains  to  be  solved.  The  analytical  background  of  this  algorithm  is 
not  new.  The  best  account  of  it  may  be  found  in  [11.  The  numerical,  microcomputer 
oriented,  algorithmic  approach  presented  in  this  report  is  new,  however,  and  highly 
efficient. 

The  other  algorithm  is  based  on  direct  optimization  of  the  LLF  associated  with  the 
PDF  (“).  It  may  be  somewhat  slower  to  execute  than  the  other  one,  but  it  has  the 
advantage  of  avoiding  the  problems  created  by  multiple  roots  of  the  equation  h(X)  =  0. 
As  a  matter  of  fact,  if  the  estimation  problem  has  a  solution  at  all,  h(X)  has  at  least 
two  distinct  zeros. 

In  many  practical  situations,  the  shift  parameter  is  known  to  have  a  specific  value 
X*.  In  this  situation  the  general  three-parameter  problem  reduces  to  a  two-parameter 
one.  It  is  easily  solved  by  means  of  the  ML  derivative  equations  approach.  The  ML 
optimization  algorithm  in  the  two-parameter  case  uses  a  LLF  which,  for  programming 
efficiency,  is  maintained  as  a  function  of  three  variables  in  an  artificial  way  (see  Sec. 
V1II.4). 

Moment  estimates  of  the  parameters  (two  or  three  depending  on  whether  the  shift 
parameter  is  known  or  not)  are  easily  obtained.  They  are  used  as  starting  values  for 
the  ML  optimization  algorithm. 


Questions  of  hypothesis  justification  and  goodness  of  fit  are  cuiside  the  framework 
of  this  report  and.  thus,  are  not  addressed  here. 

A  program  diskette  will  be  made  available  upon  request. 


II.  THE  DIFFUSION  CHARACTER  OF  THE  LOGNORMAL  DISTRIBUTION 
It  is  well  known  [2],  [3]  that  the  function 


”1  2 

■(x.t)  =  (y^b)  exp  -  (b“'x)  .  X  eR. 


b  =  b  (t)  = 


r  n  1/2 

L^TTtJ  .  r  =  0 


2o<x  (1  -  exp  -  2rt) 


r  ^  0  . 


is  the  delta  function  initial  condition  (at  x  =  0,  t  =  0)  solution  of  the  autonomous 
Fokker-Planck  equation 


[A(x)i(x.t)j^  [D(x)z(x.t)]j^-2^(x.t)  -  0  .  xeP.  .t>0 


A(x)  =  cx  >  0 
D  (x)  =  -  tr  X 


(diffusion  coefficient)  , 
(drift  coefficient)  . 


The  subscripts  in  (II. 3)  signify  partial  derivatives.  For  r  =  0,  (11. 3)  reduces  to  the 
standard  heat  equation. 

With  t  in  (II. 2)  considered  as  a  parameter  the  function  (II. 1)  becomes  the  normal 
(Gauss)  PDF  with  scale  parameter  b  . 
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The  transformation  x  =  log  'y  generates  from  (II.  1 ),  after  y  has  been  replaced 
by  X  ,  the  function 


g(x.t)  =  f  (log?f”'x.t)  bxj  exp  - b"^  log^  x  (11.4) 

with  b  =  b(t)  given  by  (11.2).  (For  details  on  this  kind  of  PDF  transformations  the  reader 
is  referred  to  [4. (5. 6)].)  For  any  fixed  t  >  0  and  with  b  =  /l  a.  log  2f  =  p,  this 
function  takes  the  form  (**)  with  X  =  0.  The  shift  parameter  X  in  (*•)  is  immaterial  in 
this  context. 

The  same  transformation  applied  to  the  differential  equation  (11.3)  leads  to  the 
autonomous  Fokker-Planck  equation 


[a“(x)  w(x.t)]j^  -  [d"(x)  w(x.t)]j^  -  w^(x,t)  =  0  , 

A^Cx)  =  (xX^  , 

D^Cx)  =  [«  -  r  log  x]  X  . 

of  which  the  function  g(x.t)  defined  in  (11.4)  is  a  solution.  As  a  matter  of  fact,  the 
function 

v"(x.t;y)  =  t)x)  exp-b~2  [ log  v"’ x  -  log  tr"' y] 


with  b  =  b(t)  again  defined  by  (11.2),  is  the  delta  function  initia’  condition  (at  x  =  y.  t  = 
0)  solution  of  (II. 5).  Relative  to  applications,  for  example  in  nucleation  evolution 
processes  [5,  Secs.  5.3.1,  5.3.31,  it  is  of  interest  to  note  that  the  scale  parameter 

p  =  log  y)  exp  -  rtj 

depends  on  time  t  unless  the  drift  parameter  x  is  zero. 
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I'l.  DATA  SETS 


In  population  estimation  problems  the  parameter  vector  P  =  (X.a.p)  of  the 
lognormal  distribution  class  (**)  is  to  be  specified  from  a  given  set  of  N  observations. 
(In  general  terms,  the  observations  represent  not  necessarily  distinct  i.i.d.  sample 
values  of  a  random  variable  which  is  assumed  to  be  distributed  according  to  some  PDF.) 

It  may  be  assumed  that  the  given  observations  form  a  set  of  h  ^  N  distinct 
elements  xi  <  X2  <  ...  <  xn  •  located  on  a  coordinate  axis  of  x  ,  each  of  absolute 

M 

frequency  (integral  multiplicity)  f  ^  1,N=  Z  f  >  1. 

^  v=l 

In  histogram  estimation  problems,  the  parameter  vector  P  =  (X.a.p)  of  (**)  is  to 
be  specified  from  a  given  sei  of  absolute  frequencies  fgy  ^  0  associated  with  adjacent 

class  intervals  of  the  form  [a+(v-l)Aa,  a  *  vAa)  (v=l . M),  a  e  R  ,  Aa  >  0,  of  equal 

length  Aa  ,  located  on  a  coordinate  axis  of  x  .  Actually,  the  vector  P  is  to  be 
specified  form  the  set  {(xy.  fa^))  of  the  coordinate  pairs  (xy,  fgy)  with  Xy  taken  as 
the  midpoint  of  the  class  interval  with  number  v  ,  i.e.,  Xy  =  a  ♦  (v-l/2)Aa.  In  this 
sense  the  histogram  estimation  problem  is  a  curve  fitting  problem. 

There  is  a  variant  of  the  histogram  estimation  problem  which  is  based  on  given 
data  sets  {(xy.  fav))  '^'th  arbitrarily  spaced  observations  X]  <  X2  ...  <  x^. 

it  is  convenient  to  have  a  uniform  notation  to  cover  all  of  these  estimation 
problems.  Therefore,  it  will  be  assumed  that  the  given  data  are  of  the  general  form  of 
a  set  of  ordered  pairs  {(xy,  fay))  (v=  1 ....  n)  with  x  j  <  X2  <  ...  <  x^  ,  and  M  s  N  = 

M 

Z  f  .  In  population  estimation  problems  risN.  fav^  1.  whereas  in  histogram 
v=l 

estimation  problems  M  <  N,  fgy  ^  0  (v=2.  ...  .  M-1).  fg]  s  1.  fgM  ^  '■  ^nd  Xy  =  a  *  (v- 
1/2)Aa. 

The  shift  parameter  X  .  in  either  case,  now  is  necessarily  restricted  to  the  half- 
interval  X  <  X  ] . 

For  the  analytical  description  of  the  estimation  problem  it  is  also  convenient  to 

-1  ^ 

introduce  the  relative  frequencies  f  =  N  '  f  .  Z  f^  =  '■ 

<3y  ^3, 


It  may  happen  that  a  value  X**  is  known  for  the  shift  parameter,  either  from 
physical  or  other  evidence.  This  case,  which  reduces  the  three-parameter  problem  to  a 
two-parameter  one  with  only  a  and  p  unknown,  is  incorporated  in  the  general 
estiication  algorithms. 
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IV.  THE  MOMENTS 


The  first  three  population  moments  for  the  PDFC**)  are  given  by 


M,  =  J  xfO<)dx  =  . 

X 

«>  2 

X 

00  3 

^3  "  (a>+2)(6)-i)^  ,  £i)  =  exp  . 

X 

M2  >  0,  M3  >  0  Since  to  >  1. 


The  corresponding  sample  moments  are,  with  the  notation  conventions  introduced 
in  Section  ill, 


v=i 


Also  of  interest  are  the  moment  coefficients  of  skewness. 


and 


-3/2 

33  ^  m3 


It  should  be  observed  that  the  number  33  may  turn  out  to  be  nonpositive  for  a  given 
data  set  whereas  ©<3  ,  which  has  been  derived  from  the  theoretical  moments,  is 
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always  positive.  Therefore,  if  s  0,  moment  estimation  of  the  lognormal  parameters 
is  impossible.  The  reader  is  eferred  to  [1.  Sec.  51  for  additional  remarks. 

Estimation  of  the  three  parameters  X,  a.  and  p  for  the  PDF  (♦*)  can  now  be 
accomplished  as  follows.  Since  o<3  is  independent  of  p  and  X  ,  the  equation  a3  = 
<X3  ,  in  which  a3  is  a  known  number,  leads  to  the  equation 

pCo)  =  =  0  . 

Application  of  Descartes'  Rule  of  Signs  [6.  Lemma  6.21  shows  that  pCo))  has  exactly  one 
simple  positive  zero  (Oq  which,  obviously,  is  greater  than  unity.  It  provides  the 
moment  estimate  for  the  shape  parameter  o  , 


1/2 


The  equation  m2  =  II2  .  then  leads  to 


(IV.  0 


Pq  =  1  {logm2-log[o>^(<^o'0]} 

Finally,  the  equation  m^  =  Mi  gives 

X  =  m ,  -  u  exp  y  . 

0  1  0  ^  ^0 

Thus,  the  vector  of  the  parameters  estimated  by  means  of  the  moments  is  Pq 
Oq,  Pq)- 


(1V.2) 


=  (Xq, 


If  the  value  X**  of  the  shift  parameter  X  is  already  known,  estimates  Oq  and 
Po  and  y  can  be  obtained  from  the  first  two  moments.  Elimination  of 

immediately  leads  to 

“0  = '  *  ^"2  ("’1  ■  ^”)  ■ 

With  this  value  of  (Oq  the  moment  estimates  for  o  and  y  follow  from  (IV. l )  and 
(IV. 2).  respectively. 
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V.  THE  LOGARITHMIC  LIKELIHOOD  FUNCTION 

In  general  terms,  but  on  the  basis  of  the  notational  conventions  adopted  in 
Section  HI.  the  likelihood  function  of  a  random  variable  X  for  which  there  exists  a 

data  set  {(Xy,  f3y))  (v=l . M)  and  which  is  assumed  to  be  distributed  according  to  a 

PDF  f(x:P).  which  depends  on  a  parameter  vector  P  =  (X.  p^.  p2  .  ...).  X  =  shift,  is 
defined  by 


L(P)  = 


f 


2  •  •• 


r^n[<v 


(see.  e.g..  [7.  Sec.  12.5].  [8.  Sec.  5.4].)  For  the  lognormal  PDF  class  (**)  with  P  =  (X.a.p) 
and  with  the  abbreviation 


log  (Xy  -  X)  =  py(X)  .  X  <  xi  . 


the  function  L(P)  takes  the  particular  form 


L(P)  =  (2Tr)‘N/2  exp-N 


109  O  ♦  Vy(X)  *  (2(?2)  ^  i^f  (pyW-p) 
v=i  v=i 


The  function 

<t.(p)  =  N'’  iog[(2Tr)'^/2^pJ 


log  a 


V=1  V=1 


(V.l) 


is  called  the  logarithmic  likelihood  function  (LLF)  of  the  lognormal  distribution  class 


The  problem  now  is  to  specify  numerical  values  for  the  parameters  X.  a.  and  p 
for  which  (**)  approximates  the  given  data  optimally.  The  maximum-likelihood  principle 
asserts  that  the  optimal  parameter  values,  if  they  exist,  are  the  coordinates  of  the 
point  P  =  (S.,  a,  jj)  in  the  open  parameter  space  P  ^  {X  <  x;  .  a  >  0.  p  e  R)  at  which 
the  LLF  <}>(P)  takes  its  global  maximum.  The  openness  of  the  space  P  is  essential  in 
this  context  as  will  be  seen  in  Section  VI. 
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VI.  riAXIhUM-LIKELlHOOD  PARAMETER  ESTIMATION 


The  first  approach  is  based  on  the  derivative  equations  which  result  from  (V.1). 
A  necessary  condition  on  the  parameters  is  that  they  satisfy  the  three  partial 
derivative  equations  8<}>/9X  =  0  .  d^/da  =  0.  and  9<}>/3|i  =  0.  They  are  of  the  form 


=  if^exp-p  -cr'^  £'v(pv''') 

OA  y-  ,  y-  ] 


exp-p^  =  0 


(VI,  1) 


f 


9<^  _ 
9a 


(VI.2) 


^  =  a'2  £  f  = 


(V1.3) 


Equations  (VI.2)  and  (VI. 3)  can  be  written  as 


-- 

v=i 


(VI.4) 


P 


i 

ij“  1 


(VI.5) 


If  the  shift  parameter  X  is  known  these  two  equations  immediately  represent 
the  solution  of  the  ML  parameter  estimation  problem  on  the  basis  of  the  derivative 
equations.  In  this  situation,  equation  (VI. l)  becomes  vacuous,  and,  for  a  given  data  set 
{(xy.  known  X  =  X**,  equation  (IV.5)  determines  ji  ,  and  (VI. •4),  with  p  =  p  , 

determines  a  [8,  Sec.  5.2.2,  Table  5. 1]. 

If  the  shift  parameter  X  is  considered  as  unknown  (and  thus,  together  with  a 
and  )i  ,  to  be  estimated  from  the  given  data),  equations  (VI. 4)  and  (VI.5)  must  be 
supplemented  by  equation  (VI.  1 ).  However,  by  means  of  (VI. 5)  and  (VI. 4),  p  and  o 
can  be  eliminated  from  (VI.  1)  so  that  only  one  equation  in  X  remains.  It  is  convenient 
to  perform  this  elimination  in  the  LLF  <j>(P)  given  in  (V.  1 )  directly.  To  simplify  the 
notation  the  following  functions  of  X  are  introduced: 
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A(X) 


P?.(X) 


>  0 


C(X)  =  ^  =  P  ■ 

V=1 

E(X)  -  -C'(X)  =  ^  f^exp-p  (X)  >  0  . 
v=l 


F(X)  =  -  ^a'(X)  =  exp-p^(X)  . 

2  v=l 

Then  (VI. 4)  takes  the  form 

cj2  =  a(a)  -  c2(x)  >  0  . 

Positivity  of  this  expression  follovv's  from  Tchebychef's  inequality  [9,  Thm.  431. 
Elimination  of  p  and  cr  (by  means  of  the  expressions  given  above)  from  <|>(X,cf,p) 
results  in  the  function 

HX)  =--j[log(A-cO  *2C*  i] 

For  9(.X)  to  take  a  maximum,  the  derivative  equation 


=  0 


with 


h(x)  -  (A  -  c2  -  c)E  F 

must  be  satisfied.  Therefore,  since  A  -  c2  >  o,  the  maximum-likelihood  estimate  X] 
for  the  shift  parameter  X  must  be  a  root  of  the  equation 

h(X)  =  (A  -  c2  -  C)E  ♦  F  =  0.  (VI. 6) 

The  following  should  now  be  observed.  Since  C(X)  i  -  oo  as  XT  X)  =  min  (xy), 
since  A(X)C"2(x) -*  f|"l  >  1  as  X  T  x|  ,  and  since  log  |C(X)(  increases  more  slowly 
than  |C(X)|  ,  the  rewritten  version 
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l[2!og|C(X)|  -  log(A(>0c'2(X)-l)  *  2Cm] 


for  $(X)  shows  that  <1>(X)  t  «  as  X  T  X].  As  a  consequence  of  this  result  it  is  seen 
that  the  openness  of  the  parameter  space  P  =  {X  <  xi  ,  a  >  0.  p  e  R  )  is  essential  for 
the  identification  of  the  location  of  the  global  maximum  of  (}>(P)  in  the  interior  of  P 
(see  the  last  paragraph  of  Section  V).  Another  consequence  is  that,  if  (I>(X)  takes  a 
maximum  at  all.  it  must,  in  the  direction  of  decreasing  X  .  take  a  minimum  first. 
Consequently,  if  the  equation  (VI. 6)  has  a  root  at  all.  it  has  at  least  two.  The  correct 
root  relative  to  the  estimation  problem  is  the  second  one  in  the  direction  of  decreasing 
X  . 


If  X]  is  the  correct  root  of  (VI. 6)  the  estimates  pi  and  c?i  for  p  and  o 
are  obtained  as 


r  T 

=  C(X,)  .  a,  =  [a(x,)-c2(x,)J 
respectively. 

The  second  ML  parameter  estimation  procedure  for  the  lognormal  distribution 
class  is  based  on  direct  optimization  in  the  parameter  space  P  of  the  LLF  <|>(P)  given 
in  (V.  1).  In  optimization  procedures  it  is  customary  to  minimize  the  objective  function. 
Therefore,  instead  of  maximizing  <t>(P),  the  function  -<j>(P)  vnll  be  minimized.  The 
resulting  coordinates  of  the  minimum  of  -<|>(P)  will  be  denoted  by  X2 . 02.  P2- 


Since  both  ML  parameter  estimation  algorithms  are  based  on  the  LLF  <|>(P)  it  is 
obvious  that  the  resulting  parameter  estimate  vectors  P]  and  P2  should  agree  within 
the  tolerances  associated  with  the  different  numerical  procedures.  Agreement  of  the 
vectors  Pj  and  P2  indicates  successful  execution  of  the  ML  estimation  process. 

In  the  numerical  implementation  of  the  optimization  ML  estimation  approach  the 
moment  estimation  vector  Pq  =  (Xq.  Oq,  Pq)  (see  Section  IV)  is  used  as  the  starting 

vector  for  the  iteration  processes.  If  the  value  X**  of  X  is  known,  the  starting 
vector  is  Pq  =  (X**,  Cq,  Pq). 
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Vll.  USER'S  GUIDE 


This  section  provides  information  on  some  of  the  numerical  aspects  of  the 
estimation  algorithms.  Information  about  their  use  is  given  on  the  diskette. 

1.  Preparation  of  Data 

According  to  Section  HI.  the  data  input  for  the  algorithms  is  a  set  {(xyi  ^av)^ 
(v=l,  ...  ,  M)  of  coordinate  pairs,  ordered  with  respect  to  the  abscissa  values  such  that 
xi  <  X2  <  ...  <  Xf-|.  Unless  the  data  are  supplied  in  this  particular  form,  a  subroutine 
does  the  ordering  and  frequency  determination  for  given  population  data.  For  histogram 
data,  a  subroutine  calculates  the  Xy's  for  given  a  e  R  and  Aa  >  0. 

For  computational  reasons,  it  is  convenient  to  translate  the  original  Xy  values 
along  the  coordinate  axis  of  x  . 

(a)  Shift  Parameter  Unknown 

The  desired  translation  is  accomplished  by  addition  of  -X]  to  each  Xy  so 
that  the  new  xj  =  0.  The  advantage  of  this  operation  is  to  have  a  universal  upper 
bound,  namely  0.  for  the  unknown  (translated)  shift  parameter  X  .  The  shape  and 
scale  parameters  a  and  p  are  not  affected.  The  algorithms  calculate  the  parameter 
values  X,  c,  and  p  and  return  the  true  parameter  values  X^  =  X  +  X],  and  pi 
relative  to  the  original  coordinate  axis. 

(b)  Shift  Parameter  Known 

Let  X**  be  the  known  value  of  the  shift  parameter  X  .  The  desired 
translation  is  accomplished  now  by  addition  of  -X**  to  each  Xy.  As  before,  o  and  p 
are  not  affected.  The  algorithms  calculate  the  values  a  and  p  with  shift  equal  to 
zero,  and  return  the  true  values 
Xt  =  X**,  0]  and  p  . 

2.  noment  Estimation 

Moment  estimation  is  straightforward  according  to  the  formulas  given  in 
Section  IV,  provided  33  >  0.  In  the  three-parameter  case,  the  equation  p(w)  =  0  is 
solved  to  find  its  root  ooq  >  1.  This  is  done  by  means  of  root  bracketing  and  a 
modified  version  of  Brent's  method  [10,  Chap.  7.31. 

3.  Derivative  Equations  ML  Estimation 

In  the  three-parameter  case  the  equation  h(X)  =  0  (VI. 6)  is  solved.  Bracketing 
of  the  desired  root  X]  is  achieved  by  calculation  of  h(Xy)  along  the  search  sequence 
Xy*  1  =  (1 .61 8)Xy  (v  =  0,  1 ,  ...).  with  Xq  = -0.05.  until  h  changes  from  negative  to 
positive  values  at  two  successive  points.  Once  X]  has  been  bracketed  the  modified 
Brent's  method  is  used  to  locate  it. 
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With  \]  calculated  (or  with  X**  known)  the  estimates  and  a\  are 
calculated  from  (VI. 5)  and  (VI. 4).  respectively. 

4.  Qotimi7ation  ML  Estimation 

Optimization  of  the  LLF  (V.l)  is  accomplished  by  means  of  Powell's  method  [10. 
Chap.  10.51  in  both  the  three-  and  two-parameter  cases.  The  vector  Pq  =  (Xq,  cTq,  Pq) 

(or  Pq  =  (X**,  Oq,  Pq))  of  moment  estimates  is  used  to  start  the  minimization 
routine  on  the  function  -<|>(P).  Once  Powell's  tolerance  has  been  achieved  the  algorithm 
continues  by  minimizing  the  function 

'{'(P)  =  (exp  10)  [^-1  ♦exp -10  (<|>(P)-<|>(P))] 


in  which  P  is  the  final  vector  from  Powell's  routine  on  -^(P).  Use  of  the  "amplifier" 
function  'l'(P)  eliminates  the  effects  of  shallowness  of  -<|>(P)  near  its  minimum  and 
increases  the  accuracy  of  the  numerical  values  of  the  coordinates  of  the  minimum  of 
-<t>(P). 


If  X  is  known  to  have  the  value  X**,  the  LLF  (|>  reduces  to  a  function  of  only 
two  variables.  To  avoid  additional  subroutines  for  this  situation,  (j>(X**,a,p)  is 
maintained  as  a  function  of  three  variables,  a,  p,  and  X,  by  addition  of  the  term 
(X-X'‘)2. 

5.  Freouencu  Calculation 

For  histogram  parameter  estimation  an  auxiliary  program  is  provided  which 
calculates  two  sets  of  expected  absolute  frequencies.  One  set  corresponds  to  the 
moment  parameter  estimates,  the  other  to  the  ML  estimates.  (It  should  be  remembered 
that,  under  successful  execution  of  the  two  ML  estimation  algorithms,  the  derivative 
equations  and  the  optimization  estimates  agree.)  In  both  calculations  the  logarithm  of 
the  PDF(**)  is  used  to  obtain  the  relative  frequencies  at  the  points  Xy  .  They  are 
subsequently  converted  into  absolute  frequencies. 

To  provide  a  basis  for  goodness-of-fit  statements,  chi-square  values  for  the 
individual  class  intervals  are  also  calculated. 
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VIII.  EXAMPLES 


To  demonstrate  the  versatility  of  the  algorithms  presented  in  this  report  five 
examples  are  offered.  They  contain  population  and  histogram  estimates  with  shift 
parameter  unknown  or  known.  Accompanying  tables  are  given  in  Section  IX. 

Example  1 

N  =  87  observations  of  annual  24-hour  maximum  rainfalls  (in  points)  at  Sidney, 
Australia,  over  the  period  1859-1945  are  given  in  Table  1.1  [111.  Moment  and  ML 
parameter  estimates  for  the  three  unknown  parameters  are  given  in  Table  1.2.  The 
subscripts  0,  1.  and  2  refer  to  moments,  derivative  equations,  and  optimization, 
respectively. 

Grouping  into  histogram  absolute  frequency  data  has  been  performed  on  the  data 
of  Table  1.1  in  five  different  ways,  6^  .  as  displayed  in  the  second  columns  of  Tables 
1.4  Gy  (v  =  1 . 5).  The  first  column  shows  the  class  interval  numbers  v  . 

Table  1.3  contains  the  class  interval  data  for  the  groupings  Gy  (v  =  1 . 5). 

The  resulting  parameter  estimates  for  the  grouped  histogram  data  are  shown  in 
Tables  1.5  Gy  (v  =  1 . 5). 


The  estimated  parameter  values  are  used  to  calculate  the  expected  absolute 
frequencies  from  the  PDFC**).  For  the  moment  estimates  they  are  given  in  the  third 
columns  of  Tables  1.4  Gy,  and  for  the  ML  estimates  in  the  fifth  columns.  Next  to  the 
calculated  frequencies  in  Tables  1.4  Gy  are  shown  the  chi-square  values.  As 
mentioned  in  the  Introduction,  no  significance  will  be  attached  to  them  within  the 
framework  of  this  report. 

Example.  2 

N  =  20  random  observations  generated  from  a  lognormal  population  with  X  = 
100,  a  =  0.4,  and  p  =  log  50  [1.  Sec.  10]  are  displayed  in  Table  2.1.  The  three 
corresponding  sets  of  estimated  parameters  are  given  in  Table  2.2. 

Since  it  is  known  how  the  data  of  Table  2.1  have  been  generated  it  is  possible 
to  process  them  under  the  assumption  that  the  value  X**  =  100  for  the  shift  parameter 
is  known.  The  resulting  two-parameter  estimates  are  shown  in  Table  2.3. 

This  example  is  on  the  diskette  as  RANDSAM.DAT. 

Example  3 

This  example  deals  with  a  total  of  N  =  37  observations  of  frost-days  in  Munich 
in  the  months  of  April  over  the  period  1930-1966  [12,  Ex.  42,  Table  31]  which  contain 
M  =  18  distinct  ones.  The  observations  together  with  their  frequencies  are  given  in 
Table  3.1. 
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The  resulting  parameter  estimates  are  shown  in  Table  3.2. 


Example  4 

A  total  of  N  =  885  observations  of  rainfall  totals  (in  inch)  for  sets  of  four 
consecutive  months  at  Camden  Square.  London,  over  the  period  1870-1943  [13.  Ex. 

7.51 1]  have  been  grouped  into  M  =  17  class  intervals  [a  +  (v-1)  Aa.  a  +  vAa)  with  a  =  2. 
Aa  =  I.  The  data  are  given  in  Table  4.1. 

The  calculated  parameter  estimates  are  contained  in  Table  4.2. 

The  expected  frequencies  calculated  from  the  parameters  of  Table  4.2  together 
with  the  corresponding  chi-square  values  are  shown  in  Table  4.1 

It  is  reasonable  to  assume  in  the  present  example  that  the  probability  of  no 
rainfall  at  the  Camden  Square  location  over  the  given  period  is  zero.  In  other  words,  it 
is  reasonable  to  assume  that  the  shift  parameter  X  is  known  and  that  X  =  x"  =  0. 

The  corresponding  two-parameter  estimates  are  given  in  Table  4.3. 

Table  4.4  shows  the  expected  frequencies  for  the  two-parameter  case. 

This  example  is  on  the  diskette  as  BROOKS.DAT. 

Example  5 

This  example  concerns  weekly  precipitation  sums  observations  (in  10"2  inch)  at 
Kwajalein.  IS.  for  the  summers  1949-1958  (14.  p.  47.  Ex.  IC).  The  total  number  of 
observations  is  N  =  130.  They  have  been  grouped  into  M  =  17  class  intervals  [a  + 
(v-l)Aa.  a  +  vAa)  with  a  =  0.  Aa  =  50.  Table  5.1  shows  the  given  data.  The  estimated 
parameter  values  are  displayed  in  Table  5.2.  The  corresponding  frequencies  are  given  in 
Table  5.1. 

The  estimates  for  a  and  p  under  the  assumption  that  the  shift  parameter  is 
zero  are  shown  in  Table  5.3.  The  expected  frequencies  in  this  case  are  displayed  in 
Table  5.4. 
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IX.  TABLES 


15/(16  Blank) 


TABLE 

1.1. 

Annual  Maximum  Twenty- 
at  Signey,  Australia, 

Four  Hour 
1859-1945 

Rainfalls 

(in  Points) 

370 

752 

662 

190 

375 

393 

395 

301 

566 

616 

445 

403 

324 

260 

690 

392 

389 

281 

489 

753 

569 

204 

425 

648 

433 

64S 

485 

466 

263 

275 

836 

S66 

295 

434 

330 

310 

236 

487 

273 

60S 

301 

423 

441 

637 

157 

477 

364 

363 

362 

571 

177 

475 

708 

441 

652 

337 

325 

342 

316 

653 

486 

414 

418 

374 

258 

188 

484 

459 

469 

239 

391 

302 

780 

322 

360 

335 

263 

216 

339 

579 

230 

333 

154 

1  105 

330 

192 

420 

TABLE  1.2.  Three-Parameter  Estimates 


LUG-NORMAL  DISTRIBUTION  CALCULAfiONS 
Data  file  is  RAINAUS.UAT 
Number  at  parameters  is  3 


•»  -ii-  -;!■  -ii-  -(HJ  is-  ->  *  if  -tf  <!•  *  *  ■»  ■»  <)■  -i:-  h;-  -:!■  -k-  -;!■  ■«  -i 


* 


THE  MOMENT  ESTIMATES  FOR  THE  LOG-NORMAL  DiSTRIGUMUN 


LambdaC  0j 
Sigma( 0j 
Mu(  0) 


-  .  722546E-*-001^ 
. 3440S7E+000 
. 615229E+001 


■!f  ■i.'- >  ->  it  a  it  it  it  *  -;i-  it  it  it  it  it  it  it  -if  it  it  it  it  it  it  a  it  it  -;)■  it 


THE  DERIVATIVE 
Lambclai  1 J  * 
Sigma(  1 )  = 

MuM  J 


ESTIMATES  FOR 
- . 394497E+001 
. 400045Et000 
. 5984b2E+001 


THE  LOG-NURMAL  DIST 


^  ^  ^  ^  JL  ^  ^  ^  Ji.  ^  JC  Jt  Jt  Ot.  jL  -L  1  V  V 


THE  PUWELL  ESTIMATES  tUH  THE  LOG-NURMAL  D 1 h T 
LambdaC 2J  -  - . 3e9397E+00 1 

Sigma(2j  =  .400101E+000 

Mul  2)  =  .  59d437E<-001 


it  -c-  «•  it  -M-  -:)•  it  it  it  h;-  it  a  it  -  t  it  it  it  it  it 
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TABLE  1.3. 


Class  Intervals  and  Numbers  of  Classes  for  Groupings  G„ 
(v»l,  ...  ,5) 


<5, 

^3 

^4 

a 

150 

125 

150 

100 

125 

Aa 

100 

100 

50 

75 

75 

n 

10 

10 

20 

14 

14 

TABLE  1.4G]^.  Grouped  Frequency  Data  for  Grouping  Gj 


LOG-NOHMAL  DISTHIBUTIUN  CALCULATlONb 
Date  file  is  Ingdsetl.dat 


Moment  Estimate  Max  Likelihood  tst 


V 

f(  abs) 

f  (  ca 

1) 

X«*2 

f  (  c  a 

L) 

x-;^<s-2 

1 

1  1 

1  0 

.  64 

.01 

1  0  , 

.  83 

.  00 

2 

23 

21 

.  05 

.  16 

24  , 

,  62 

.  1 1 

3 

23 

20 

.  97 

.20 

20  . 

,  95 

.  20 

4 

10 

i4 

.  92 

1 . 63 

13. 

.  38 

.  65 

b 

10 

8 

.  88 

.  14 

7  . 

,72 

.  68 

6 

4 

4 

.  77 

.  12 

4  , 

,  30 

.  02 

7 

4 

2 

.42 

1 . 02 

2  . 

,  38 

1  .  1  1 

8 

1 

1 

.  19 

.  03 

1  , 

.  32 

.  08 

9 

0 

.  50 

.58 

,  75 

.  75 

10 

1 

.  28 

1 . 07 

.  43 

.  77 

Totals 

07 

05 

.71 

86 

.68 
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TABLE  1.4G2.  Grouped  Frequency  Data  for  Grouping  G2 


LOG-NORMAL  DISTRIBUTION  CALCULATIONS 
Data  file  is  Ingd8et2.dat 


Moment 

Lstimate 

Max  Likelihood 

L  s  t 

V 

f (  abs) 

f ( cal) 

fCcal) 

2 

1 

e 

7 . 56 

.  03 

7.10 

.  1  1 

2 

10 

18.90 

.  04 

20.67 

.  35 

3 

25 

2  1.69 

.  50 

22 . 23 

.  35 

4 

1b 

16.76 

.  03 

16.09 

.  00 

5 

7 

10.41 

1.11 

9.70 

.  75 

6 

7 

5 . 68 

.  30 

5 . 34 

.51 

7 

3 

2 . 88 

.01 

2.82 

.01 

8 

2 

1 . 39 

.  27 

1 . 46 

.  20 

9 

0 

.66 

.  66 

.75 

.75 

10 

1 

.31 

1 .58 

.  39 

.  97 

Total  6 

87 

86.23 

86.55 

TABLE  I.4G3. 

Grouped 

Frequency  Data 

for  Grouping  G3 

LOG-NORMAL  DISTRIBUTION  CALCULATIONS 
Data  file  is  lngdset3.dat 


Moment 

Estimate 

Max  Likelihood  Est 

V 

f(  abs) 

f(  cal) 

X*»2 

f  (  ca  i  ) 

X-*!-<»2 

1 

6 

3.62 

1 .56 

3.29 

2.23 

2 

5 

7.13 

.  64 

7 . 54 

.  85 

3 

8 

10.00 

.  40 

10.78 

.  72 

4 

15 

11.33 

1.19 

11.94 

.  7b 

5 

12 

11.13 

.  07 

11.36 

.  04 

6 

1  1 

9 . 92 

.  12 

9.62 

.  14 

7 

10 

8 . 25 

.  37 

7 . 96 

.51 

8 

0 

6 . 53 

6.53 

6 . 22 

6 . 22 

9 

5 

4 . 99 

.  00 

4.71 

.  02 

10 

5 

3.70 

.45 

3 . 50 

.  64 

1  1 

3 

2 . 70 

.  03 

2.57 

.  07 

12 

1 

1  .  94 

.  4b 

1  .  87 

.  40 

13 

3 

1 . 38 

1.91 

1 . 35 

2.01 

14 

1 

.  97 

.  00 

.  98 

.  00 

15 

1 

.68 

.  15 

.  70 

.  13 

16 

0 

.48 

.  48 

.51 

.51 

17 

0 

.  33 

.  33 

.  37 

.  37 

18 

0 

.  23 

.  23 

.  26 

.  26 

19 

0 

.  16 

.  16 

.  19 

.  19 

20 

1 

.  1  1 

6 . 97 

.  14 

5 . 32 

Totals 

87 

85.61 

86.06 

TABLE  I.4G4.  Grouped  Frequency  Data  for  Grouping  G4 


LOG-NORMAL  DISTRIBUTION  CALCULATIONS 
Data  file  is  lngdset4.dat 


V 

f(  absj 

Moment 
f(  cal) 

Est imate 

X*">2 

Max  Likelihood  Est 
f(calj  X-!:-^-;-2 

1 

2 

2.51 

.  10 

1 . 88 

.  0  1 

2 

9 

8 . 90 

.  00 

9.19 

.  00 

3 

15 

14.69 

.  00 

15.85 

.  05 

4 

20 

16.54 

.72 

17.03 

.  52 

5 

13 

14.42 

.  14 

14.28 

.  1  1 

b 

8 

10.81 

.  73 

10.41 

.  56 

7 

7 

7 . 34 

.  02 

6 . 98 

.  00 

8 

6 

4 . 66 

.  38 

4 . 45 

.  54 

y 

3 

2 . 84 

.01 

2 . 75 

.  02 

10 

2 

1  .  68 

.  06 

1  .  67 

.  06 

1  1 

1 

.  97 

.  00 

1  .  00 

.  00 

12 

0 

.  56 

.  56 

.  60 

.  60 

1  3 

0 

.  32 

.  32 

.  36 

.  3b 

14 

1 

.  18 

3.74 

.  22 

2 . 84 

Totals 

87 

86 . 62 

8b  .  65 

TABLE  I.4G5.  Grouped  Frequency  Data  for  Grouping  G5 


LUG-NOHMAL  DISTRIBUTION  CALCULATIONS 
Data  file  is  lngdset5.dat 


Mome  n  t 

Est ima  t  e 

Max  Likelihood  Est 

\j 

f ( ab  9 ) 

fC  cal) 

X"  -  2 

f  C  cal) 

X  2 

1 

6 

4 . 37 

60 

4.18 

.  79 

2 

6 

12.35 

1 . 53 

12.53 

1  .  64 

3 

20 

16.99 

.  54 

17.26 

.  44 

4 

17 

16.31 

.  03 

16.40 

.  02 

5 

16 

12.83 

.  78 

12.77 

.  82 

6 

4 

8 . 97 

2 . 76 

8 . 87 

2 . 68 

7 

6 

5 . 84 

.  00 

5 . 76 

.  0  1 

8 

4 

3 . 64 

.  04 

3 . 59 

.  05 

9 

3 

2.21 

.  29 

2.19 

.  30 

10 

1 

1  .  32 

.  08 

1.31 

.  07 

1  1 

1 

.  78 

.  06 

.  78 

.  06 

12 

0 

.  46 

.  46 

.  47 

.  47 

1  3 

0 

.  27 

.  27 

.  28 

.  28 

1  4 

1 

.  1  6 

4 . 43 

.  1  7 

4.18 

Totals 

87 

86 . 49 

86.56 
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TABLE  1.5Gi.  Three-Parameter  Estimates  for  Grouping  Gj 

LCG-NORMAL  DISTRIBUTIUN  CALCULAT lONb 
Data  file  is  LNGUSt.  T  1  .  DAT 
Number  of  parameters  is  3 


■tt- -it  •iHi- -ii- ■'.i- -M- Hr -> •»  ■i;  -ij  •{(•  h;- hs- -> -t- » « * h} h:- Hr-  -if 


THE  MOMENT 
Lambdal 0)  = 
SigmaC  0) 
MuC0) 


ESTIMATES  FOR  THE 
- . 8389 1 0E+002 
.342081E+000 
.  61a104E-^001 


LOG-NORMAL  DISTRIBUriON 


*  Hi-  *  H}  Hi-  Hi  Hi-  Hr-  Hi-  Hi-  Hi-  Hi-  Hi-  Hi-  Hi-  Hi-  Hi  Hi-  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  -U-  Hi  Hi  Hi  Hi  Hi  Hi  Hi  H 


THE  DERIVATIVE 
L  ambda(  1 J  = 

S igmat  1 ) 

MuC  1 J 


ESTIMATES  FOR 
. 731777E+002 
.505160E-^000 
. 574925E+001 


(HE  lug-normal  Ulbl 


^  -ii-  -i'c  "iJ-  -.r  **}•  -IJ-  -i'c  iz  ie  iz  iz  i'i  i'z  iz  iz  w  -.'z  •/:  z'z  zz  *. 


THE  POWELL 
Lambda( 2J  = 
&igma( 2) 
Vu(2) 


ESTIMATES  FOR  THE 
.731776E+002 
.  505  160E-^000 
.  574925E-t-001 


LOG-NORMAL  DIST 


Hi  Hi  Hi  Hi  -'i  Hi  Hi  Hi  Hi  Hi  Hi-  -ii  Hi  Hi  HT-  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hr  Hi  Hi  Hi  Hi  r  Hi  -.i  ,i  -.i  Hi  H:  Hi 


TABLE  1.5G2.  Three-Parameter  Estimates  for  Grouping  G2 


LOG-NORMAL  DISTRIBUTION  CALCULATIONS 
Data  file  is  LNGDSET2.DAT 
Number  of  parameters  is  3 


j'.  j/.  ^  jj,  .y.  .‘j,  /j.  .y.  .‘j,  ^  >*.  .;l  .it  ••  -ij.  .jt  .y. 


THE  MOMENT 
Lambda( 0)  = 

S igma(.  0)  = 

Mu(  0) 


ESTIMATES  FOh  THE 
-  .  122396E-t-003 
.  31764, HE  +  000 
.  62599bE  *-00  1 


L  0  G  —  N  Ij  H  M  A  i.  U  i  b  F  H  1  d  U  i  i  1 N 


J/-  -JJ.  -‘J.  -j;-  Jj.  Jy.  J.  -jj.  zf,  >t  ^  ^  -•<.  iw  vt.  jj.  .«t  .y.  Jv  ,*•  w;  jh.  j-.  .  ^  .  ..  .  .  .  h. 


THE  DERIVATIVE 
Lambda!.  1  }  = 

SigmaC  1 J  = 

Mu(  1  J 


ESTIMATES  FOR 
-.141  9b4E-H102 
.  398  1  59L-i000 
. 60 1 29 1 £+00 1 


the  LOG-NURMaL  DIaT 


Hi  Hi  Hi  Hr  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi 


’t  Ji.  Zi.  Z',  J 


THE  POWELL  ESTIMATES  FOR  THE  LOG-NORMAL  Ulb 
LambdaC 2)  =  - . 1 4 1 964E +002 

SigmaC  2J  •=  .3901S9b+000 

MuC  2}  =  .601291E  +  001 


Hi  -;i  Hr  Hi  -,i  Hi  -:i  -  r  Hi  Hr  Hi  Hr  Hi  Hi  Hi  Hi  Hi  H.-  Hr  Hi  vi  Hi  Hi  Hi  Hi  Hr  Hi  Hi  Hi  Hi  Hi  Hr  Hi  H.  Hr  Hr  Hr  Hr  :i  H;  Hr  H.  H.  H,  Hi  i  Hi 
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TABLE  I.5G3.  Three-Parameter  Estimates  for  Grouping  G3 

LOG-NORMAL  OISTRIBUTION  CALCULATIONS 
Data  file  is  LNGDSET3.DAT 
Number  of  parameters  is  3 


^  .'j.  ^  ^  ^  JL  ^  ^  .jc  ^  ^  ^ 


THE  MOMENT  ESTIMATES  FOR  THE  LOG-NURMAL  DISTRIBUTION 


Lambdal  0j 
Sigma( 0J 
M  u  i  0 ) 


-  .  574054E+00iJ 
. 3S9288E+000 
.81 1734L+001 


J*.  JJ.  JJ.  ^  Jj.  J*.  ^  Jt  .ij.  rfj.  ^  rffc  JU  J 


THE  DERIVATIVE 
Lambdat  1 )  = 

bigmai.  1 J 
Mu(  1) 


ESTIMATES  FOR 
.  1b8447E  +  002 
. 422920E+000 
.592ei0E+001 


THL  LOU-NOHMAL  UlS 


Wj.  .'J.  ./t  JJ.  Jj.  .JJ.  JJ.  JJ.  -‘t  .>C  ^  Jl.  .*4  ^  ^  -"4  ^  >4  .'4  ^  -J4  J.  Ji-  J4  V  _  J'.  i'.  .  4  J  .  .  .  J'.  .  J-  .J  V 

THE  POWELL  ESTIMATES  FOR  THE  LOG-NORMAL  DTST 
Lambdat2)  »  .158448E+002 

Sigma( 2J  =  .422920E+000 

Mu(2j  =  .S92810E+001 

-:)■  Hi’  Hi-  •;(•  •«•  li-  Hi-  ->  •»  Hi  •»  -ii-  ■M-’M-  *  •!}  «-  ■«■  •»  -if-  -K-  •>!•  -S!-  Hi-  «•  -K-  ■>}•  Hi-  Hi  •>:•  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi 


TABLE  I.5G4.  Three-Parameter  Estimates  for  Grouping  G4 


LGG-NORMAL  distribution  CALCULATlDNb 
Data  file  is  LNGUSET4.DAT 
Number  of  parameters  is  3 


JJ  JJ.  Jj  vj.  .;<•  J4  J4  JH,  .*4  JJ.  J4  .-J  J4  JJ.  JJ,  J4  JJ  JJ,  JJ  JJ  .^J  ,^J  JJ  JJ  ■/.  ^4  Jj  JJ  .;4  jJ  Jj  .y.  Jj  v  .  .  •.  J4  .  .  .  v  . 

THE  MOMENT  ESTIMATES  FOR  THE  LOG-NURMAL  DibrRiBUTiuN 

Lambda( 0J  =  - . 1 0473 1 E+003 

Sigmal,  0.)  =  .  326079EH-000 

Mu(.0)  -  .622973E+001 


JJ  JJ  JJ  jj  JJ  jj  JJ  .JJ  ^4;-  -ij  ^J  JJ  -JJ  ^  ^  -iJ-  -{J  ^  -iV  -ir  ■  v  ■%<  \  -  %'  *”*  -  -  - 


THE  DERIVATIVE  ESTIMATES  FOR 
Lambda( 1)  “  - . 265338E+002 
Sigma(  1J  =  .  382665EH-000 
Mu(1)  =  .605208E+001 


THE  LOG-NURMAL  DIST 


H;-  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  Hi  ..  ■.  ,i 

THE  POWELL  ESTIMATES  FOR  THE  LUG-NUHMal  IJibT 


LambdaT  2 J  = 
S igmat  2 )  = 
Mu(2J 


-  .  265343E’+’002 
. 382665E+000 
.  605209EH-00  1 


■  "Jr  iy  O’  iy  ■>r  O'  i'r  -yr  >>  -ii-  •>«*  ir 
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TABLE  I.5G5.  Three'Parameter  Estimates  for  Grouping  G5 


LOG-NORMAL  DISTRIOUTION  CALCULATIONS 
Data  file  la  LNGOSt TS . OAT 
Number  of  parameters  is  3 

^  JC  -V  J'.  JV  -‘V  Jj.  J;.  ^  J}.  J.  .'t.  i-.  -JJ.  J'.  .  .‘j.  >j.  .‘v  .y.  J'.  _‘t  •  .>4.  -  ,  •  -><-  J;  J-.  ;•  ....  -• 

THE  MOMENT  ESTIMATES  f OR  THE  LOG-NURMAL  DiblRIBUliUN 
LambdaC 0j  =  - . 436032E+002 

Sigma(0j  =  .370006E+000 

MuT 0)  =  .607b71E+001 


^  ^  Jt  JJ.  ^  ^  -  Jt  ^  JJ.  >1  -‘J.  ^  jL  .fC  ^  jy.  j<.  .;t.  jj,  jj.  ^  Jj.  -;L  .;j.  J 

THE  DERIVATIVE  ESTIMAIES  FOR  THE  LOG-NORMAL  DISi 

LambdaC  1j  =•  -  .  246  1  30E  +  002 

SigmaCI}  =  .384999E+000 

Mu( 1)  «  .603142E+001 

^L  JJ.  JL  JC  JL  -^4  VL  -JV  Jt  JJ.  .jC  «*J.  Jt  V.  M.  JJ.  JU  ^  .  i',  If.  Jf.  .  •.  Jf.  .•'.  .  .»{.  ^ 

THE  HOWELL  ESTIMATES  FOR  THE  LOG-NORMAl  UIST 

LambdaC  2)  =  -  .  246 1 62L +  002 

SigmaC 2)  *  .384997E+000 

MuC2)  =  ,603143E  +  001 


%;■  -i;- ■!)•  -M- -i;- -:s- ■» a- •«■  •«■  •«•  -i;- <•  <!■ « a- a- ii-  :> c-  i:-  a <- -;i- 


23 


TABLE  2.1.  A  Random  Lognormal  Sample 


146.290 

133.143 

132.971 

164.304 

145.788 


184.101 

144.328 

156.680 

128.709 

155.369 


135 . 880 
166.475 
174 . 600 
153.070 
201.415 


TABLE  2.2.  Three-Parameter  Estimates 


LQG-NOHMAL  DlSTHiBUTlON  CALCULATIONS 
Data  file  is  fiANOSAM.DAT 
Number  of  parameters  is  3 


■ii-  *  4!-  it  -i)-  ■>  -ii-  -it  it  -a-  ■>  it  it  it  it  it  it  it  it  it  it  <■  .> 


IHL  MOMENT 
LambdaC  0J 
Sigmal,  0J 
Mut  0J 


tSTlMATES  KOH  THE 
. 727727t+002 
.241002fc+000 
. 434708t+00 1 


LUG- NUHMAL 


UiST 


j'.  JL  .y.  .y.  jy.  jy.  jy.  .jj,  .{t  ty.  ^y.  Jj.  jL  J*.  Jy.  -;w  .y.  ^ 


IHE  DERIVATIVE 
LambdaC  1 J  = 

Si gmai  1  j  = 
Mull) 


ESTIMATES  LOR 
.  1  1772  IE +  003 
. 60431 3E+000 
.  337317E  +  001 


IHE  LOG-NUHMAL  L 


v;.  jy.  jy.  jy.  .‘t  .;j.  .;j.  >yf  .JC  ^  Jt  jy.  .. 


JV  J‘.  -•>  .JV  J\ 


•-  J'_  .K.  JJ.  J*-  .'t 


IHE  POWELL 
Lambda(  2 ) 
SigmaC  2} 
Mu(  2) 


ESTIMATES  PUR  IHE 
. 1 17784E+003 
. 605e9yE+000 
. 337039E+001 


LOG-NOPIMAL 


Dis  r 


^  -Yf  "Jr  ■ir  -ii-  -[•  m'  ^  "ii*  ■ii-  "Vr  •i^  -i^  ■vJ"  ^  -ir  -iJ-  <{■  ^}-  *i5-  ^}-  -it  *{’•  vi-  -ii-  -ii-  -Yr  •vi-  ■ 


2  4 


127.211 
1 37 . 338 
131 . 375 
168.554 
157.238 


i  b  U  I  i  U  N 


iST 


TABLE  2.3.  Two-Parameter  Estimates 


LOG-NORMAL  DISTRIBUTION  CALCULATIONS 
Data  file  is  RANOSAM.DAT 
Number  of  parameters  is  2 


*  » -;f  -if  a-  -if  -i;-  if  -if  if  if  if  if  if  if  if  if  if  if 


-  .>•-  >-  ->  -'<•  -  - 


’’HE  MOMENT  ESTIMATES  (-QR  THE  LOG-NORMAL  UISIHIHUIION 


Lambda!.  0.) 
Sigma(.  0j 
Mut.  0) 


.  1  00000E  +  003 
.  127181E  +  000 
.  501978E  +  001 


if  if  :f  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  if  -.f  if  if  if  i.-  if 


THE  DERIVATIVE 
Lambda(.  1 )  * 

Sigma! 1 )  = 

Mu!  1J 


ESTIMATES  EUR 
.  1000008  +  003 
. 385203E+000 
. 389015E+001 


THE  LUG  ’NURM/^l  01  S  ! 


,‘J.  J J.  JJ.  J'.  .JJ.  -J'.  J/.  JL  ^  .-L  Jj.  J'.  J'.  . 


.  i'.  J'-  < 


THE  HOWELL 
Lambda(  2) 
Sigma( 21 
Mu(2l 


ESTIMATES  E OH  THE 
=  .100000E+003 
=  .366203E+000 
=  .389015E+001 


log-normal  UIST 


if  if  if  if  if  if  if  if  if  if  if  if  if  i;- if  ♦  if  if  «  if  if  if  if  if  if  if  if  if  if  if  if  if  if 


Jj.  .jj,  .fj.  .jL  JL  -;j.  j  . 
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Q  n  >4  Q 


TABLE  3.1.  Number  of  Frost-Days  in  the  Months  of  April,  1930-1966 
and  Their  Frequencies 


1 

3 

4 
6 

10 


2 

3 

B 

8 

10 


2  2  2 

3  4  4 

S  S  6 

0  0  9 

12  17 


2 

4 

6 

9 


TABLE  3.2.  Three-Parameter  Estimates 


LOG-NORMAL  0  I S TR i BO f I ON  OALCULATIONS 
Data  file  is  LNDSET4.DAT 
Number  of  parameters  is  3 


,{4.  jj.  ^4.  ^4.  ^  ^  •iHi*  •ii'  i'f  -f'-  •Ji*  -i'c  “JJ  5“  %> 

the  moment  ESTIMATES  FOR  THE  LOG-NORMAL  DiSTRlBUTlON 
Lambda(0j  =  - . 639374E+00 1 

Bi£mat0)  =  .288822E+000 

Mu( 0)  =  .24444BE+001 


-if -;nfj 


.4  . 


}.  Jt  J,  J'-  ^4.  .-fc  J 


4  ^  .'4.  •  ‘j*  •  •  •' 


THE  DERIVAIIVE 
L  ambda(  1 J  = 

8 i gma( 1 )  = 

Mut  1  } 


ESTIMATES  FOR 
- . 30208bE+001 
. 407262E+000 
. 207493E+001 


THE  LOG-NORMAL  OIS 


J4.  Jt  .  .  j’. 


i'c  v*  ^  if 


-;4.  .‘v  .K.  ^  ^  .;j.  JL  . 


THE  POWELL  ESTlMATLb  TOR  I  HE 


L  ambda( 2  j 
Sigmat  2J 
MuC  2J 


- . 302085E+00 1 
. 407262E+000 
. 207493E+00 1 


LUG-NuHMmL 


U  ib  T 


J'.  JL  J'.  iJJ  J4.  JJ.  H.  ^  >t.  -;t.  JJ.  t/.  JJ.  .•«.  Jj.  .;4.  .Jt  J;-  j  .  .T- 
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C\J  O  O) 


TABLE  4.1.  Rainfall  Totals  (in  Inch)  for  Sets  of  Four  Consecutive 
Months,  Camden  Square,  London,  1870-1943 


LOG-NORMAL  DISTRIBUTION  CALCULATIONS 
Data  file  is  brooKs.dat 


Moment 

Estimate 

Max  Likelihood  bst 

V 

f(  abs) 

f(  cal) 

X*»2 

f(  cal) 

X*->2 

1 

9 

4 . 77 

3.75 

4 . 69 

3.97 

2 

15 

19.92 

1 .21 

19.70 

1.16 

3 

48 

50.43 

.  12 

50 . 35 

.  1  1 

4 

81 

09.11 

.74 

09.19 

.75 

5 

142 

120.66 

3.77 

120.05 

3.70 

6 

125 

133.54 

.55 

133.73 

.  57 

9 

129 

126.52 

.05 

126.63 

.  04 

e 

109 

106.15 

.08 

106. 16 

.  08 

9 

76 

80.89 

.  30 

60 . 84 

.  29 

10 

56 

57.07 

.  02 

57.01 

.02 

1  1 

36 

37.86 

.09 

37.80 

.  09 

12 

22 

23.88 

.  15 

23.84 

.  14 

13 

13 

14.46 

.  15 

14.44 

.  14 

14 

12 

0 . 47 

1 . 47 

0.46 

1  .40 

15 

5 

4.83 

.01 

4.83 

.01 

16 

4 

2 . 69 

.63 

2.69 

.63 

17 

3 

1 .40 

1  .57 

1 . 48 

1  .  57 

Totals 

80S 

002.73 

802.75 

27 


TABLE  4.2.  Three-Parameter  Estimates 


LOG-NORMAL  DISTRIBUTION  CALCULATIONS 
Uata  file  is  BROOKS. DAI 
Number  oT  parameters  is  3 


Jt  jj.  ^  .jj.  jj.  .'j.  Jj.  .  j.  ^  .jt  jj.  jt  j'.  jt  .y-  .jj.  -y.  -  .  j.  y.  -  j. 

THE  MOMENT  ESTIMATES  FOR  THE  LOG-NORMAL  UlSTRIBUliUN 
LambdalBJ  =  -  .  4388  1  6E  +  00  1 

SigmaC0J  =  .215054E+000 

Mu( 0)  =  .2B3126E+001 


.‘t  Jt  .‘V  jt.  jj.  jj.  .‘j.  jt  -*j.  >j.  .IL  Jt.  jj.  ^  > 


THE  DERIVATIVE 
Lambda( 1 J  = 
Sign.a(  1 ) 

Mu(  1 J 


ESTIMATES  FOR 
- . 431653E+001 
. 2 1 6064E+000 
. 2b2544E+001 


THE  log-normal  UIsT 


•)i-  ->  J.}-  *!f-  4;- 


i'r  i'y  \c  ■{:-  v-  •«-  -{i*  ir  •{<•  i'r  if  •>  i'r  i'f  i(  if  if  /f  if  ir  if  i 


i'.  i'r  •: 


THE  POWELL 
Lambda!.  2.1  = 

Sigma( 2j  = 
Mu(2j 


ESTIMATES  FOR  THE 
- .431654E+001 
.216064E+000 
.2S2644E+001 


LOG-NORMAL  DIST 


i;-  •«•  -if  ■i,'-  -it-  i;-  ->  -i;-  <;■  -e;-  <>■  •«■ 


«  « •«•  a  -a-  ■«• 


/L 
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TABLE  4.3.  Two-Parameter  Estimates 


» 


% 


LOG-NGHMAL  OISTRIBUTJGN  CALGULATIDNS 
Data  file  is  BROOKS. DAT 
Number  of  parameters  is  2 


<!•  -it-  ■*>  ■!!■  a  a  v-  *  -ii-  -m-  ■JS'  •"  -*  ■"  ■i*'  "■  '*■  '■ 

IHt  MOMLNT  LSTIMATBS  FOB  THt  LOG-NOHMAL  UlblHlBUliON 
Lambda!  0J  =  .0O0000t+000 

Sigma(0J  =  .321?0bt+000 

Mut 0J  =  .20e539E+001 


.-t  .'c  -  j.  .'C  Ji.  Jt  jj.  ^  ^  ^  jc  .‘j.  J-.  >j.  ^  j'.  j}.  -/♦ 


THt  UEHI\/AriVE 
Lambda!  1 J  = 
Sigma(  1 J  = 

Mu(  1J 


ESTIMATES  FOR 
.  000000E-*-000 
.  341S4SE  +  000 
. 200 1 32E+00 1 


THE  log-normal  UIST 


X.  Jv  Jj.  -JJ*  -{J  -fr  -i;- 


JC  Ji,  .'i.  Ji.  J*  • 


THE  F^OWELL 
Lambda( 2) 
Sigma( 2) 

MuC  2) 


ESTIMATES  FOB  THE  LOG-NORMAL 
.000000E+000 
. 341S45E+000 
.20ai32t+001 


J'. 

DiST 


.;!■  i;.  -a-  •«••«•  *  -K-  -» ■:!•  -i:-  •»  •'>  -t-  *  ii  «  <!■  ^  ii  i-  •-  ii-  "  "  "  - 
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TABLE  4.4.  Expected  Frequencies  for  Two-Parameter  Estimates 


LOG-NOHMAL  DISTRIBUTION  CALCULATIONS 
Data  file  is  brooks.dat 


Moment 

Estimate 

Max  Likelihood  Lst 

V 

f (  abs) 

f(  cal) 

X-H-H-2 

f ( cal) 

X*-if2 

1 

9 

.58 

122.01 

1 . 23 

49 . 09 

2 

15 

10.88 

1 .56 

15.57 

.  02 

3 

40 

47.45 

.01 

55 . 07 

.91 

4 

01 

99.01 

3.20 

102.35 

4 . 45 

5 

142 

135.60 

.  30 

131.76 

.  80 

6 

12B 

143. 13 

2.30 

135.25 

.  78 

7 

129 

127.51 

.02 

119.83 

.  70 

8 

109 

101.29 

.59 

96.14 

1 . 72 

9 

76 

74.31 

.  04 

72  .  £2 

.  22 

10 

56 

51.54 

.  39 

51.41 

.41 

1  1 

36 

34 . 35 

.  08 

35 . 47 

.  0  1 

12 

22 

22 . 24 

.  00 

23.08 

.  15 

1  3 

13 

14.12 

.09 

15.81 

.  50 

14 

12 

0.03 

1.14 

10.34 

.  27 

IS 

5 

5.47 

.  04 

6.71 

.  43 

16 

4 

3.37 

.  12 

4 . 33 

.  02 

17 

3 

2.06 

.43 

2.79 

.  02 

Totals 

885 

861.73 

879.95 
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TABLE  5.1.  Weekly  Precipitation  Sums  (in  10“^  inch),  Kwajalein,  IS, 
1949-1958 


LOG-NORMAL  OISTRIBUTION  CALCULATIONS 
Data  file  is  fq.Kwa 


Moment 

t  s  t ima t  e 

Max  Likelihood  tst 

V 

f  (.  ab s) 

f( cal) 

X««2 

f(  cal) 

X  “  "‘■2 

1 

8 

10.20 

.48 

8 . 84 

.  08 

2 

25 

15.68 

5.54 

22.26 

.  34 

3 

21 

18.00 

.47 

23.01 

.  18 

4 

15 

17.50 

.  30 

18.98 

.  84 

5 

18 

15.36 

.45 

14.40 

.  06 

6 

9 

12.51 

.  99 

10.73 

.  28 

7 

8 

9.73 

.31 

7 . 89 

.  00 

0 

6 

7.33 

.24 

5.01 

.01 

9 

5 

5.40 

.  03 

4 . 30 

.  12 

10 

4 

3.92 

.  00 

3.20 

.  20 

1  1 

3 

2.82 

.01 

2.40 

.  15 

12 

1 

2.02 

.51 

1  .  82 

.  37 

1  3 

1 

1  .43 

.  13 

1  .  39 

.  1  1 

14 

1 

1 . 02 

.  00 

1 . 07 

.  00 

15 

2 

.72 

2.25 

.  83 

1 . 64 

16 

2 

.  52 

4.28 

.65 

2 . 00 

17 

1 

.37 

1 .09 

.51 

.46 

Totals 

130 

124.60 

128.19 
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TABLE  5.2.  Three-Parameter  Estimatep 


LOG-NORMAl  UI&THiBUriGN  GALCULATIONG 
Data  file  is  PQ.KWA 
Number  of  parameters  is  3 

j;.  .  J.  .‘J.  .y.  JJ.  -JJ.  .y.  ,»j.  jj.  .y.  Jy.  jj.  Jy.  .}.  .jj.  .*1.  jj,  -‘i,  -y.  j-  v-i. 

THE  MOMENT  ESTIMATES  FOR  THE  LOG-NORMAL  DISTRIBUTION 
Lambda!  D)  =  -  .  1  BWASyL +  0k3 J 

bigma( 01  =  .410B79E+000 

MulDl  =  .593b04L+001 

.  .  J-.  ,  .  -JJ.  jJ.  aj.  .-L  JJ,  Jj.  JJ  Jt  JJ.  Jj.  J^  wj.  J^.  JJ.  JL  -  J  J^  JJ.  JJ.  JJ.  Jj.  Jj.  .  J.  JJ.  .j.  .  j.  J  ..  y.  ^y.  jy.  jy.  .y.  .  .  _  .  v  ..  .  . 

IRE  UEHIVAnvE  ESIIMATEb  f  UR  IRE  LUG-NUHMAL  Ulbl 
Lambda!  1)  =  -.2169L'1E-'002 

bignia!  1.)  =  .  6942b0E  •-000 

Mu!  1 J  =  .53006 1 E  +  00 1 

THE  PQWElL  ESTIMATES  FOR  THE  LOG-NORMAL  DIET 
Lambda! 2J  =  - . 2 1 692 i E +002 

Sigma! 2]  =  .6y4260E+000 

Mu!  2J  =  .530061E  +  001 


Jt  Jt  .JJ.  JJ.  JJ.  Jt  JJ.  JJ.  jfc  Jv  J,J.  ^  JJ.  JV  jJ.  ^  ^j.  -jt  J^.  ^  Jj.  Jt  jJ.  JJ.  JJ  Jt  JJ  .  J  jj  JJ.  Jv  J^-  JJ  Jj.  JJ  JJ  jj  JJ  JJ  j j 
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TABLE  5.3.  Two-Parameter  Estima-es 


LUU-NOHMmL  L)iS  I  HIBU  r  IUN  BALBIJLA  r  IQNB 
Data  tile  is  FQ.KWA 
Number  o+'  parameters  is  2 

Jj.  Jj.  .  JJ.  .;j.  j;.  .;}.  J  .  j'.  jj.  jj.  .jj.  jj.  .jj.  Jj.  -*V  .jv  JJ.  .jt  .jj.  Jj.  .•«.  - 1  .‘j.  /•.  -«•.  vC  .»(.  .jt  Jj-  J/.  -■■.  ;j.  . 

THt  MOMENT  ESllMAPES  FOR  THE  LOG-NORMAL  DlblfllBUFION 
Lambda(.0J  =  .000000E  +  000 

bigma(,0J  =  .b79SJ2E  +  000 

Mu(0J  =  .521ie9E+00? 

.  J*.  Jv  J;.  j-..  .  v  Jj.  .Jt  -^4  .'j.  Jj.  .‘j.  .JJ.  -y.  .jj.  ^  .y,  jj.  jj.  .y,  .•».  j".  ./  /-  .  .  .  .•».  jj.  vv  -  - 

IHE  OERiVATIVE  ESTIMATES  FOR  THE  LOG-NORMAL  UIST 
Lambda!,  1J  =  .000000L+000 

bi^mal.  1)  =  .ai9Vb3£+000 

Mui  1 J  =  . 5 14332L  +  00 1 

J  .  jy.  .;j.  -;j.  4*.  .  i.  .*l  .*j.  .H  jj.  jj.  .jv  .-t  .jt,  .h  .;j.  jv  ;j.  .‘L  .  .  /j.  -;t  jy.  .• 

THE  POWFILL  E&riMAiES  FOR  7  HE  LOG-NORMAL  U  i  t- T 
Lambda!,  2j  =  .000000E+000 

bigma!.2j  -  .8iy7b3E+000 

Mu(2J  =  .514332E+001 

J/.  j;.  j;.  .  l  rft.  wj.  vj.  .)t  JL  v'.  y.  vc  .M  .«<.  Jv  ^  .y.  .jv  jy»  .y.  /c  .y.  j.  jy.  Jy.  . y  .JJ.  .J;.  jt. 
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TABLE  5.4.  Expected  Frequencies  for  Two-Parameter  Estimates 


LOG-NORMAL  DISTRIBUTION  CALCULATIONS 
Data  file  is  fq.Kwa 


Moment  tstimate  Max  Likelihood  Lst 


V 

f ( ab  s ) 

f ( cal J 

X*«-2 

f ( cal) 

1 

8 

2 . 02 

17.66 

8 . 04 

.  00 

2 

25 

21.35 

.  62 

25 . 39 

.  0  1 

3 

21 

26 . 08 

.  99 

23.51 

.  27 

4 

15 

2  1.02 

2.13 

18.07 

.  52 

5 

18 

16.25 

.  19 

13.30 

1  .  66 

6 

9 

11.64 

.  60 

9 . 74 

.  0b 

7 

8 

8 . 25 

.  0  1 

7.17 

.  1  0 

8 

6 

5 . 85 

.  00 

5 . 34 

.  08 

y 

5 

4.17 

.  1  6 

4 . 03 

.  24 

1  0 

4 

3.01 

.  33 

3 . 07 

.  28 

1  1 

3 

2.19 

.  30 

2 . 37 

.  17 

12 

1 

1.61 

.  23 

1  .  85 

.  39 

13 

1 

1.19 

.  03 

1 . 45 

.  14 

14 

1 

.  89 

.  0  1 

1.16 

.02 

15 

2 

.  67 

2.60 

.  93 

1 . 24 

1  6 

2 

.51 

4.29 

.75 

2 . 09 

17 

1 

.40 

.  92 

.  6  1 

.  25 

Total  3 

130 

127.91 

126.77 
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